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ABSTRACT 

The recently observed X-ray synchrotron emission from four supernova 
remnants (SNR) has strengthened the evidence that cosmic ray electrons are 
accelerated in SNR. We show, that if this is indeed the case, the local electron 
spectrum will be strongly time-dependent, at least above roughly 30 GeV. The 
time dependence stems from the Poisson fluctuations in the number of SNR 
within a certain volume and within a certain time interval. As far as cosmic ray 
electrons are concerned, the Galaxy looks like actively bubbling swiss cheese 
rather than a steady, homogeneously filled system. 

Our finding has important consequences for studies of the Galactic diffuse 
gamma-ray emission, for which a strong excess over model predictions above 
1 GeV has been reported recently. While these models were relying on an 
electron injection spectrum with index 2.4 - chosen to fit the local electron flux 
up to 1 TeV - we show that an electron injection index of around 2.0 would a) 
be consistent with the expected Poisson fluctuations in the locally observable 
electron spectrum and b) explain the above mentioned gamma-ray excess above 1 
GeV. An electron injection index around 2 would also correspond to the average 
radio synchrotron spectrum of individual SNR. We use a three-dimensional 
propagation code to calculate the spectra of electrons throughout the Galaxy 
and show that the longitude and latitude distribution of the leptonic gamma-ray 
production above 1 GeV is in accord with the respective distributions for the 
gamma-ray excess. 

We finally point out that our model implies a strong systematic uncertainty in 
the determination of the spectrum of the extragalactic gamma-ray background. 



Subject headings: gamma rays: theory - cosmic rays - ISM: supernova remnants 
- acceleration of particles 
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1. Introduction 



As was first observed by OSO-3 ( [Krausiiaar et al. 1972| ), the dominant feature of 



tlie liigli-energy 7-ray sky is tlie intense emission from the Galactic plane. Later the 
complete SAS-2 ([Fichtel et al. 1975 ) and COS-B ( Mayer-Hasselwander et al. 1982 ) data 



gave evidence for a correlation between the 7-ray emission and the spatial structures of 
the Galaxy. The intensity distribution and the spectral form of the emission have led 
to the consensus that the diffuse 7-ray radiation is primarily produced by interactions 
between Galactic cosmic ray particles and the interstellar medium, and to a small extent by 
unresolved Galactic point sources ( [Bloemen 1989| , [Strong 1995| ). The EGRET observations 



of the Magellanic Clouds have shown that cosmic ray nucleons in the energy range below 
100 GeV are almost certainly galactic ( [Sreekumar et al. 1993| ), while the observations made 



with the OSSE and COMPTEL instruments aboard CGRO have provided strong evidence 
that cosmic ray electrons are galactic ( [Schlickeiser et al. 1997 , see also i'azio et al. 1966). 



Therefore the diffuse Galactic 7-ray emission tells us about the propagation of cosmic 
rays from their sources to the interaction regions and thus complements the direct particle 
measurements by balloon and satellite experiments. 

The greater sensitivity and spatial and energy resolution of EGRET compared 
with SAS-2 and COS-B permit a much more detailed analysis of the diffuse Galactic 
7-ray emission than was possible with the earlier experiments. The spatial and spectral 
distribution of the diffuse emission within 10° of the Galactic plane have recently been 
compared with a model calculation of this emission which is based on realistic interstellar 
matter and photon distributions and dynamical balance ([Hunter et al. 199^ ), i.e. cosmic 
rays having the same spectrum and composition everywhere in the Galaxy and having an 
intensity which follows the surface density of thermal gas convolved with a Gaussian with 
dispersion a = 1.7kpc — 2.0 kpc ( [Bertsch et al. 1993|) . The distribution of the total 



intensity above 100 MeV agrees surprisingly well with the model predictions. However, at 
higher energies above 1 GeV the model systematically underpredicts the 7-ray intensity. 
If the model is scaled up by a factor 1.6, the model prediction and the observed intensity 
above 1 GeV agree well. Thus the model displays a deficit of ~38% of the total observed 
emission which depends, if at all, only weakly on location. At energies above 1 GeV around 
90% of the model intensity is due to 7r°-decay, i.e. hadronic processes, and only 10% is due 
to interactions of electrons. 

There are a number of possible explanations for this deficit: 

- A miscalibration of EGRET could cause an overestimation of the intensity above 1 
GeV. This possibility is highly unlikely. Point sources generally show power-law spectra 
without spectral hardening above 1 GeV. It would require an extreme level of cosmic 
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conspiracy for a calibration error to mimic a general power-law behavior in the spectra of 
cosmic 7-ray sources. 

- The kinematics of tt" production may be poorly understood. Detailed Monte-Carlo 
calculations have shown ( [Mori 19971) , that models based on the current knowledge of particle 
interactions do not give very different results for the 7r° spectra than do simple isobar plus 
scaling descriptions ( Permer 1986| ). It is unlikely that the cosmic ray nucleon spectrum in 



the solar vicinity is softer than that elsewhere in the Galaxy. The local cosmic ray spectrum 
samples sources within a few kpc in distance and a few times 10'' years in time. Since the 
observed deficit appears to be independent of Galactic longitude, the sources of cosmic rays 
within a few kpc from the sun would have to be different from those in the inner Galaxy 
and those in the outer Galaxy. We have also tested and verified that the uncertainties in 
the local interstellar cosmic ray spectrum below a few GeV are by far not sufficient to 
account for the deficit. The uncertainties arising from our limited knowledge of the cosmic 
ray nucleon spectrum and the nucleon-nucleon interaction kinematics can be estimated to 
be on the order of a few percent. 

- There may be unresolved point sources which contribute strongly at higher 7-ray 
energies. The only known class of objects with appropriate spectra is pulsars. Based on the 
properties of the six identified 7-ray pulsars it has been found that unresolved pulsars would 
indeed contribute mainly between 1 GeV and 10 GeV ( Pohl et al. 1997a| ). However, to 
account for all the deficit it is required that more than 30 pulsars be detectable by EGRET 
as point sources. This can be compared with less than 10 unidentified 7-ray sources which 
are not variable ( [McLaughlin et al. 1996| ) and which show pulsar-like spectra ( [Merck et 



al. 1996[) . Also, the latitude distribution of the 7-ray emission from unresolved pulsars is 



inconsistent with that of the observed emission. Unresolved pulsars will contribute 6-10% of 
the observed 7-ray intensity above 1 GeV and around 3% in the energy band between 100 
MeV and 1 GeV ( [Pohl et al. 1997a[ ), and thus they can account only for a small fraction of 
the high-energy 7-ray deficit. 

All in all the effects described above can account only for a small fraction of the deficit, 
or can add only small systematic uncertainties. In this paper we will investigate whether 
the remaining deficit of 30-35% may be caused by inverse Compton emission of high-energy 
electrons. Leptonic processes contribute only around 10% of the intensity above 1 GeV 
in the model of Hunter et al. (1997), which corresponds to ~6% of the total observed 
intensity. Thus the leptonic contribution would have to be increased to 35-40% of the total 
observed emission to explain all the deficit. 



The Galactic distribution of cosmic ray electrons is intimately linked to that of 
their sources. The recently-found evidence of X-ray synchrotron radiation from the four 
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supernova remnants SN1006 ([Koyama et al. 1995| ), RX J1713. 7-3946 ( Koyama et al. 1997| ) 



IC443 ( [Keohane et al. 1997| ), and Cas A ( [Allen et al. 1997| ) supports the hypothesis that 



Galactic cosmic ray electrons are accelerated predominantly in SNR. X-ray synchrotron 
radiation implies TeV 7-ray emission from the comptonization of the microwave background 
at a flux level which depends only on the average magnetic field strength (|Pohl 1996|) , and 



indeed the detection of the remnant SN1006 at TeV energies has been announced recently 
( P?animori et al. 19971 ). Interestingly, there is no clear observational proof that the nuclear 
component of cosmic rays is accelerated likewise in SNR. The acceleration of cosmic ray 
nucleons in SNR should lead to observable flux levels at TeV energies ( Prury, Aharonian 



and Volk 1994|) , however with a spectrum different from that of the leptonic emission. 



The generally tight upper limits for TeV emission from the nearest SNR ( |Lessard et al. 



1995| , [Buckley et al. 1998| ) are in conflict with simple shock acceleration models for cosmic 



ray nucleons in SNR. 

Most of the radio synchrotron spectra of SNR can be well represented by power-laws 



with indices around a ~ 0.5, corresponding to electron injection indices of s ~ 2.0 ([Green 



1995 ). This is in accord with predictions based on models of particle acceleration ( [BlandforH 



and Eichler 19871) . However, it is different from the electron injection spectral index of 
s = 2.4 which has been inferred from the locally observed electron spectrum ( [Skibo 1993[ ) 
and which subsequently has been used in the model of Galactic 7-ray emission of Hunter 
et al. (1997). The contribution of cosmic ray electrons to the Galactic 7-ray spectrum at 
high energies depends strongly on their injection spectral index. If the acceleration cut-off 
energy is high enough, the leptonic 7-ray emission may even dominate at TeV-PeV energies 
( Porter and Protheroe 1997[) . A change in the electron injection index by 6s = 0.4 could 
increase the inverse Compton emissivities at a few GeV by an order of magnitude or more. 

Let us suppose that for some reason the local cosmic ray electron spectrum is different 
from the average electron spectrum in the Galaxy. Then the following scenario appears 
viable: the bulk of cosmic ray electrons is accelerated in SNR with an injection index 
around s ~ 2.0. The leptonic 7-ray emission at a few GeV would be much stronger than 
in the Hunter et al. model, and it may explain a substantial fraction of the discrepancy 
between their model and the observed spectra. So if we would find a mechanism or an effect 
which would cause the local electron spectrum to be different from the Galactic average, we 
may in a second step reassess the 7-ray spectra produced by Galactic cosmic rays without 
having to assume an electron injection index of s = 2.4. 

It has been noted before that the spatial distribution of cosmic ray sources affects the 
locally observable spectra ( [Cowsik and Lee 1979| , [Lerche and Schlickeiser 1982a|) . As far 



as electron acceleration in SNR is concerned, there is no evidence that the star formation 
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activity and thus the SNR production rate in the solar vicinity is significantly less than the 
Galactic average. In the next section we will show that for cosmic ray electrons, unlike 
nucleons, the local spectra above a certain energy can deviate from the Galactic average, 
even if the spatial distribution of cosmic ray accelerating SNR is homogeneous. This is a 
result of the discrete nature of SNR both in space and in time. We will use this finding 
in the third section to model the high-energy 7-ray excess as result of inverse Compton 
emission of cosmic ray electrons, albeit with a harder injection spectrum than conventionally 
assumed. 



2. The time-dependence of the local electron spectrum 

The spectrum of cosmic ray electrons in terms of the contributions of discrete sources 
like supernova remnants (SNR) has been discussed before ( |Gowsik and Lee 1979| ). These 



authors have investigated the case of continuously active sources and have concluded that 
one needs sources situated within a few hundred parsecs of the solar system, in order 
that the energy losses of electrons do not induce a cutoff in the energy spectrum. Since 
the required number of active sources exceeds that of supernova remnants by an order of 
magnitude, SNR were found unlikely to be the only source of cosmic ray electrons between 
1 GeV and 1 TeV. In that paper the diffusion coefficient had been assumed independent of 
energy. With the usual energy dependence D oc ii;(0-3-0-6) gQj-^g Qf ^j^g statements in the 
Cowsik and Lee paper would have to be relaxed. 

In this section we will also consider the finite lifetime of SNR or other possible cosmic 
ray accelerators together with the random distribution of SNR in space and time. The 
latter induces a time dependence in the local electron spectrum at higher energies, which 
stems from the Poisson fiuctuations in the number of SNR within a certain distance and 
within a certain time interval. As we will see, the discreteness of sources does not simply 
cause a cutoff in the electron spectrum, but makes it time variable and thus unpredictable 
beyond a certain energy. 

Since effects of the discreteness of sources show up only at higher particle energy, we 
may describe the propagation of cosmic ray electrons at energies of 1 GeV to 1 TeV by a 
simplified transport equation 

(97V (9 

— - ^{bE' N)-DE^ V'N = Q (1) 

where we consider continous energy losses by synchrotron radiation and inverse Compton 
scattering, an energy-dependent diffusion coefficient DE"^, and a source term Q. The 
Green's function for this problem can be found in the literature ( pinzburg and Syrovatskii 
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1964). 

G(r, r', t, t', E') = ^ '-^^ exp - \,JL (2) 

where 

X{E, E') = T-^— (e'"^ - E"'-^) (3) 



6(1 - a 

In case of discrete sources the injection term Q is a sum over all sources. For an individual 
source showing up at time to and injecting for a time period r we can write 

= go E'-^ Q{t' - to) e(to + r - t') (4) 

Without loss of generality we can set t = and obtain 

e(f-to)e(to + r-tO exp(-g-) 
{AnKf^ {l + hEt 

where 



iV = go / dt' (5) 



A = ;S + (6) 

and r is the distance between source and observer. is the contribution to the local electron 
spectrum provided by a single source (SNR) at distance r which is (was) injecting electrons 
for a time period r starting at to- The local spectrum of electrons can now be obtained by 
summing the contributions from all individual sources. In our case the distribution of SNR 
is, for ease of exposition and computation, assumed to be a homogeneous disk of radius 
Vs = 15kpc and half-thickness Zg = 0.14kpc. Other choices for the spatial distribution 
of SNR do not impose serious changes in the results, as long as the distribution is not 
structured on sub-kpc scales. 

The numerical procedure is as follows: for each volume element 6V the expected 
number of SNR injecting cosmic rays within the time interval 6t is 

(T) 

>' tot ^m] 

where Vtot is the total volume of the source distribution and ti^j = 55 years is the inverse of 
the supernova rate ( |Cappelaro et al. 1997|) . 6t has to be the maximum lookback time in our 



model which is the sum of the SNR lifetime r and the electron energy loss time ti = ^^-^^^^^ 
at 1 GeV, the lowest energy considered here. At a given time the actual number of SNR in 
that volume element is a Poissonian random number with mean A^;^, and each of the SNR 
has a birth date which is a random number uniformly distributed within [— r — ti,0]. The 
final electron spectrum is then derived by summing over the contributions of the individual 
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SNR per volume element and summing over all relevant volume elements. We calculate 400 
such 'random' spectra and thus derive the distribution of possible spectra and their spread. 

In Fig.l we show the resulting range of local electron spectra compared with the 
observed spectra ( [Ferrando et al. 1996| , [Golden et al. 1994| , [Golden et al. 198^ , P?ang 



[1984[ , [Taira et al. 1993| ). The energy density of the ambient photon fields plus that of the 
perpendicular component of the magnetic field strength is in total taken to be 3.5 eV/cm^. 
The changes in the Compton cross section in the Klein-Nishima regime of optical and 
near-infrared photon fields are neglected. The diffusion coefficient is D = 4 ■ 10^'' cm^ sec~^ 
at 1 GeV and increases with energy to the power a=0.6, in accord with results for a 
two-dimensial diffusion model fit to the local spectra of 13 primary and secondary cosmic 
ray nuclei at rigidities between 1 GV and 10^ GV ( [Webber, Lee and Gupta 1992[ ). At these 
energies electrons and nuclei will scatter off the same turbulence, except for the helicity, 
and thus their diffusion behaviour may be expected to be similar. The injection spectral 
index of electrons is s=2.0, which corresponds to the mean synchrotron spectral index of 
individual SNR ([Green 1995[ ). 



The solar modulation of electrons at a few GeV energy has a strong effect on 
the observed spectrum. Note the difference between the spectrum observed when the 
modulation level was high ( polden et al. 1994|) and the spectrum observed during the 
passage of Ulysses over the solar south pole ([Ferrando et al. 1996[) . We have crudely 
approximated the effect of solar modulation using the force-field approach ([Gleeson and[ 
[Axford 1968|) with $ = 400 MV. Thus we probably overestimate the modulation in case of 
the Ulysses data and underestimate in case of the Golden et al. data. 

The spectra at a given time are not necessarily smooth. There is no preference for the 
usual broken power-laws or power-laws with exponential cut-offs. In fact the individual 
spectra are bumpy above ~ 50 GeV and some display step-like features. As an example we 
show a particular spectrum as dash-dotted line in Fig.l, which is slightly on the low side 
between 10 GeV and 100 GeV, where it suffers further softening before it abruptly hardens 
at 300 GeV. 

Below 10 GeV the local electron spectrum is well determined. Between 10 GeV and 
100 GeV it varies with time by a factor of 2 or 3, and above 100 GeV is completely 
unpredictable. Changes in the absolute numbers for the diffusion coefficient and the 
radiative energy losses do not change the basic behaviour, but can shift the transition 
between weak and strong variability to lower or higher energies. If the energy dependence 
of the diffusion coefficient is weaker, i.e. a < 0.6, the transition between weak and strong 
variability will be faster, and vice versa (a slower transition for higher powers than a=0.6). 
For comparison we have indicated the result for an energy dependence of the diffusion 
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coefficient, D oc E 



0.33 



As shown in Fig.l, the high-energy data for the local electron flux are in accord 
with an injection spectral index of s=2.0, though in a model with steady injection and 
a smooth source distribution these data would require an injection index of around 2.4 
( ^kibo 1993|) . Concerning the distribution of high-energy electrons, the Galaxy would look 



like Swiss cheese, with holes and regions of higher density. In the line-of-sight integrals, 
which are relevant for comparison with the EGRET 7-ray data, the averaging over holes 
and high-density regions will give the same result as a model with steady injection, however 
with source index of 2.0 instead of 2.4. At higher latitudes the line-of-sight will be so short 
that regions of low or high electron density will be resolved. The leptonic 7-ray spectra in 
direction of the Galactic poles should be relatively soft since the line-of-sight integral of the 
7-ray emissivity will be dominated by the soft local spectrum. 

The absolute electron flux is reproduced if each SNR provides an energy input of 
]^q48 gj^gg ^Yie form of electrons, which is one thousandth of the canonical value of 
Xq5i gj^gg ^Qj. j^Y^Q kinetic energy input per supernova. Taken over a lifetime of 10^ years the 
corresponding power of 5 ■ 10^^ ergs/sec is less than the X-ray luminosity of SN1006 alone. 

The time variability of the high energy cosmic ray spectrum will not be related or 
even be synchronous to variability in the flux of low-energy cosmic ray nucleons, which can 
be traced by cosmogenic instable isotopes in sediments (|Sonnett et al. 1987| , [McHargue^ 



Damon and Donahue 1995| , [Kocharov 1996|) or meteorites ([Bonino 1996| ). 



A few notes should be added. We have taken supernova explosions to be completely 
independent of each other. One may expect some level of correlation in OB associations 
and SNOBs, which would make the basic effect of time dependence even more dramatic, 
since the OB associations and SNOBs would act as single sources with longer lifetime, but 
much smaller frequency of occurrence. 

Another important point is that we have assumed that all electron sources produce 
the same spectrum. In reality this doesn't need to be the case. Some SNR will produce 
electrons with harder spectrum, and another group of SNR will provide softer spectra. 
It may be that the spectral form depends on the age of an SNR. In fact the radio data 
show that SNR do have different synchrotron spectra ( preen 1995|) . If we take the electron 



injection spectral index of an individual SNR not as a fixed number, but a random variable 
following some probability function, the time-averaged spectrum (the dotted line in Fig.Q) 
will get a positive curvature. The level of time variability on the other hand will get larger. 
The shaded regions in Fig.l, in which the spectrum is contained for 68%, respectively 95%, 
of the time, extend beyond those for fixed injection index. 
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A final note concerns secondary positrons and electrons. These particles are generated 
subsequent to interactions of cosmic ray nucleons with ambient gas, and hence the effect 
discussed does not apply and the local spectrum of secondary electrons will not vary. 
Thus the observed positron fraction will also exhibit variability anticorrelated with that 
of the primary electron spectrum. If we are indeed living in a hole in the distribution of 
high-energy electrons then the positron fraction above, say, 20 GeV will be above the level 
expected in steady injection models, if the gas density within ~ 1 kpc from the sun is not 
also sub-average. This may explain the observed positron fraction in that energy range, 
which is indeed slightly above the model predictions (Parwick et al. 1997|) . 

We have seen that the discreteness of sources of cosmic ray electrons causes a 
strong variability of the local electron spectrum at higher energies. Therefore the high 
energy electron spectrum does not prescribe our choice of electron injection spectrum in 
propagation models. 

If we consider 7-ray emission in the Galactic plane, the line-of-sight integral of the 
emissivity will correspond to an averaging over the different variability states, and hence the 
7-ray intensity calculated with time-dependent models will not differ significantly from the 
results of steady-state models. The latter are much easier to compute and thus preferrable, 
but they will in general not be able to reproduce the local high energy electron spectrum 
correctly. On the other hand, in the absence of reliable data on the position and age of all 
nearby SNR it is also impossible to calculate the time-dependent local electron spectrum 
precisely, we can only infer the level of variability. Therefore we feel that, if the acceleration 
of electrons occurs predominantly in SNR, the 7-ray emission in the Galactic plane can be 
sufficiently well described with steady-state models, provided the model is not required 
to fit the local electron data above ~ 30 GeV. In the next section we will discuss a 
steady-state model for the diffuse leptonic 7-ray emission for an injection index of 2.0. 



3. The propagation of electrons in the steady state case 

There is a richness of literature on the topic of electron propagation in the Galaxy, 
including analytical solutions for the one- or two-dimensional diffusion and diffusion- 
convection problem ( [Ginzburg and Syrovatskii 1964| , [Berkey and Shen 1969| , [Lerche a"iid 



Schlickeiser 1981| , [Lerche and Schlickeiser 1982 b| , ^ohl and Schlickeiser 19^ ). These 



solutions can be well described in their basic behavior by the concept of the catchment 
sphere ([Webster 1970|) . The energy losses prevent electron propagation farther from their 



source than a critical distance p which is defined by equality of the time scales for transport 
and for energy losses, so that the spatial dependence of the Green's function of the problem 
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is basically a function which is a constant for distances less than p and which is zero beyond 
p. If the transport is governed by diffusion and if the energy loss terms do not strongly 
depend on location, then p will not strongly depend on direction, and thus a source of 
cosmic rays at the position (x', y', z') would fill a sphere of radius p with cosmic rays. Hence 
we may separate the spatial problem and the energy problem, and approximate the solution 
to the spatial problem by a Gaussian function for the catchment sphere. The Gaussian 
function is exact at higher energies where radiative energy losses dominate, but is a crude 
approximation at very low energies, where ionization and Coulomb losses are important. 

We will include escape as a catastrophic loss term, which limits p to some maximum 
value. This approximates the effect of a sudden increase of the diffusion coefficient at a 
certain height above the Galactic plane ( [Lerche and Schlickeiser 1982b| ). We regard this a 



better description than a finite boundary with density and density gradient set to zero at 
Lhaio^ a few kpc above the Galactic plane, since for a physical escape solution the density 
outside the diffusion region relates to that inside as 



0.6 



VGevj 



and thus it is definitely not zero. With a diffusion coefficient D = DqE"-, E being the 
kinetic energy, we have 



p{E) = ^4DoE-^ (9) 

where 

^e// ~ '^esc + '^loss (10) 

We can write the differential number density of cosmic rays at position (x, y, z) coming 
from a source at (x', y', z') as 

5N ^ 



X du Q{u) exp ( - / ^ ^ , ^, I exp ( — -\ (11) 



{^pf\E\ 

' dv \ I r"^^ 

V ^ tU^)\e{^)\) ''''^ V 

where 

E is the energy loss term and Q{E) is the source spectrum at position (x', y', z'). Physically 
5iV is a propagator and it can be treated like a Green's function. Given the spatial 
distribution of sources g(a;', ?/', 2') we obtain the cosmic ray spectrum at any position (x, y, z) 

as 

N{E) = [ dx'dy'dz' 5Nq{x\y\z') (13) 
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Our method thus enables us to calculate the three-dimensional distribution of electrons 
resulting from an arbitrary three-dimensional distribution of sources. 

The individual SNR may accelerate electrons with slightly different spectra. This 
would result in a positive curvature of the composite injection spectrum ([Brecher and] 



Burbidge 1972| ). To demonstrate the effect of a possible dispersion of the injection spectral 
index in the electron sources, we assume that the injection indices for individual SNR follow 
a normal distribution 

1 / (s-so)'^ 
= — exp ; 

at the energy Eq. Radio spectral index measurements at a few GHz, corresponding to 
i?o — 5 GeV, indicate fj,s ^ 0.2 ( Preen 1995| ). Then the source spectrum of primary electrons 



Pis) = exp {-'-^^\ (14) 



IS 



1 / TP \ 0.5 ui In 

Qe = qe{m.c^y-'E-^(§-) (15) 



-^0 

For /is = this reduces to the conventionally assumed single power-law. 

The energy losses due to ionization and Coulomb interactions, bremsstrahlung and 
adiabatic cooling, and synchrotron and inverse Compton emission are well described by 



-E = 7.2- IQ-^'"^ uht 



1] E ^ eUmag ( E 



1 + -3T^ ^ + 



r714meC^ uht V 3727 mgC^ 
where we have used the following abbreviations 

Urad 



(16) 



0.75 + 



Umag 



1 ( divv \ Up 



71 

T=l + 1.54— (17) 

uh is the neutral gas density, the density of ionized gas, v the bulk velocity of electrons, 
and U the energy density of the magnetic field and the ambient radiation field in eVcm~^. 

We assume the magnetic field strength to be constant over the total volume of the 
Galaxy with i? = 10 /xG. As we shall see later, this value leads to synchrotron emission 
consistent with observations. The interstellar radiation field can be calculated from the 
respective emissivities for optical and infrared emission, and the microwave background 
emission ( [Youssefi and Strong 1991] and recent updates [Strong 1997|) . The distribution of 



ionized gas has been modelled on the basis of pulsar data ( [laylor and Cordes 1993D . The 



derivation of the distribution of neutral matter will be discussed below. For the propagation 



- 12 - 



calculation all parameters of the interstellar medium are averaged over a scale of 1 kpc to 
mimic the average environmental conditions of a cosmic ray electron during its life time. 
A Galactic wind is assumed to operate in the Galactic halo, such that adiabatic cooling 
outside the disk provides similar energy losses as bremsstrahlung inside the disk, i.e. the 
energy loss terms can be written independent of the spatial location within a catchment 
sphere. Note that the energy loss terms for neighboring catchment spheres may be different, 
since they are averages over different volumes. The radial extent of the Galaxy is taken to 
be 16.5 kpc. Note that the computer time consumption scales as radial extent squared. The 
calculation of the bremsstrahlung and inverse Compton emissivities is described elsewhere 
( IPohl 19941 ). 



The energy loss time scale, which determines the radius of the catchment sphere p, can 
be understood as an ensemble average 

_^SJdE'Q{E') J^^duii-' 



J^dE'QiE') 

where Q{E) is the electron injection spectrum. This average age can be approximated by 



Tic 



l-E 

/ du ii-^ , Ec = 2.718 E (19) 



to better than 30% accuracy except for the lowest energies. For E^ < 714r?7~^ m-eC^, 
corresponding io E ^ 100 MeV, the average age f is larger than the e-folding energy loss 
time scale. Here the energy losses will also depend explicitely on p since ioniziation and 
Coulomb interactions occur only in the gas disk. For ease of computation we will assume 
Tioss to be constant at these low energies, as if ionization and Coulomb interactions would 
not occur. This means we overestimate p but underestimate E. The two effects work in 
opposite directions but may not balance each other. We want to keep in mind that we are 
using a crude approximation at low electron energies, which for this paper however will 
have an impact only on the bremsstrahlung spectra at ^ 50 MeV. 

We neglect secondary electrons in our model. The locally observed fraction of secondary 
electrons is of order 10% ( [Barwick et al. 1997| ), but the fraction may be strongly dependent 



on position and on the propagation behavior of particles ( ^chlickeiser 1982| ). There are 
basically two arguments which may allow us to neglect the secondaries. Firstly, because of 
the energy dependence of the cosmic ray secondary-to-primary ratios, secondary electrons 
will have a production spectrum which is softer than the production spectrum of cosmic 
ray nucleons by 5s ~ 0.6, at least above 1 GeV electron energy. Thus secondary electrons 
will have a softer spectrum than primary electrons, so that their contribution to the diffuse 
Galactic 7-ray emission can not lead to a hardening of the spectrum, irrespective of the 



- 13 - 



flux. The second argument concerns the luminosity. Since the production cross sections 
for charged and for neutral pions are of the same order, and the electrons take only about 
two thirds of the pion energy, the source power supplied to secondary electrons is linked to 
the hadronic 7-ray luminosity. Only a fraction of the source power is channeled back into 
7-ray emission, since synchrotron radiation takes away some energy. Due to the kinematical 
low-energy cut-off at ~100 MeV in the secondary production spectrum and the decrease 
of hadronic 7-ray luminosity with energy above 1 GeV, the secondary contribution to 
leptonic 7-ray emission above 100 MeV will always be limited to less than 10% of the 7-ray 
luminosity due to vr^-decay ( Pohl 199^) and thus be negligible in this energy range. 



3.1. The distribution of gas 

The three-dimensional distribution of thermal material in the Galaxy has been 
determined by deconvolution of HJ surveys for atomic gas and CO surveys for molecular gas 
with the rotation curve of [Clemens (1985)|. The HI surveys include the Leiden- Greenbank 



survey ( [Burton and Liszt 1983| , [Burton 1985[ ), the Weaver- Williams survey ( [Weaver and 



Williams 1973[ ), the Maryland-Parkes survey ( [Kerr et al. 1986[ ), the high-latitude Parkes 



survey ( [Cleary, Heiles and Haslam 1979[ ), and the Heiles-Habing survey ( [Heiles and Habin^ 



1974 ). The CO data are taken from the Columbia survey ( Pame et al. 1987 ) as updated 



( pigel and Dame 1997[) . All these data are publicly available from either the ADC or the 



CDS data bank. It should be noted that none of these surveys is stray-light corrected. The 
level of uncertainty in the rotation curve, the position of the sun, and non-circular motions 
of gas is high, so that our deconvolution should be taken as a model rather than as a datum. 

The general procedere in the deconvolution process is as follows: The rotation curve 
is used to calculate the relation between distance and line-of-sight velocity, which is then 
transformed into a probability distribution for distance on the basis of the actual velocity 
resolution of the surveys and a turbulent velocity dispersion of 10 km/sec for CO and 
25 km/sec for HI. This approach tends to smear out the gas distribution along the 
line-of-sight, but relaxes most of the forbidden velocity problem. The near-far ambiguity 
towards the inner galaxy is resolved by dividing the intensity according to the amplitudes 
of Gaussian probability functions for the distribution of gas normal to the Galactic plane. 
These Gaussian probability functions are the convolutions of the local gas distribution 
functions and the spatial resolution function of the particular survey. The effective scale 
heights of gas on the far side of the Galaxy are thus systematically larger than the local 
values. As a result the gas tends to be more evenly distributed over the Galaxy than in 
other derivations ( [Hunter et al. 1997[) and the deconvolved gas distribution on the far side 
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of the Galaxy will be slightly smeared out normal to the Galactic plane. This has no impact 
on the calculation of the 7-ray emission since the column density of gas is always preserved, 
and it has also no impact on the cosmic ray propagation since our algorithm is based on the 
gas surface density, which is also preserved. 

The HJ data are scaled under the assumption of a constant spin temperature of 
Tg = 125 K. The obvious absorption features in the direction of the Galactic Center 
have been replaced by a linear interpolation between the neighboring velocity bins. The 
distribution of HJ normal to the Galactic plane is assumed to be a Gaussian of dispersion 
Zc = 0.12 + 0.023(r — 9.5) G(r — 9.5) kpc, where r is the galactocentric radius and 6 is a 
Heavyside function, with an offset according to the warping of the HJ disk ([Burton 19761). 



The CO-data are scaled with an X-factor of 1.25, which is the mean of the best fit 
values in published papers on EGRET data analysis ([Hunter et al. 1994] , Pigel et al.| 



T995i pigel et al. T996| , [Strong and Mattox 1996| , [Hunter et al. 1997[) . In the inner kpc of 



the Galaxy the X-factor is reduced to 25% of its nominal value to account for the higher 
excitation temperature and different metallicity ( [Sodroski et al. 1995[ , [Arimoto et al.[ 



[1996[ ). The noise in the CO spectra is preserved in the deconvolution process to keep the 
line-of-sight integral of the density unchanged. The vertical CO distribution is assumed 
to be a Gaussian of dispersion Zc = 0.074 + 0.03(r — 9.5) 0(r — 9.5) kpc, where r is the 
galactocentric radius and G is a Heavyside function (Pame et al. 1987| ). The position of 
the sun is assumed to be located 8.5 kpc from the Galactic Center and 15 pc above the 
plane ( Pammersley et al. 1995[) . 



For a region of ~ 20° towards the Galactic Center and towards the anticenter 
the kinematical resolution is insufficient, and the data have been edited by hand. The 
distribution of gas in these two regions is basically an interpolation between the results 
for the adjacent regions. In the anticenter region any excess over this interpolation has 
been evenly distributed over all distances, while in the Galactic Center region any excess is 
attributed to the Galactic Center. The distribution of gas in the Galactic plane is shown in 
Fig.|. 



3.2. The spatial distribution of sources 



The true distribution of SNR in the Galaxy is not well known as a result of selection 
effects and the absence of a proper distance measure to the remnants. In many papers 
( ^tecker and Jones 1977[ , Pogiel and Uryson 1988| , pioemen et al. 1993[) the cosmic ray 
distribution in the Galaxy has been estimated on the basis of a functional form for the 
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SNR distribution which fits the data for 116 remnants ( |Kodaira 1974 ). One of the general 



findings in these studies is that the overall cosmic ray distribution is too steep to explain the 
gradual slope of the 7-ray emissivity over the Galactic radius ( [Strong et al. 195^ , [Strong 
and IVlattox 1996|). 



Here we use a revised functional form for the SNR distribution ([Case and Bhattacharya 




1996|) , which fits the data for 194 remnants 

1.69±0.22 / s 

exp ( -(3.33 ± 0.37) ^ j (20) 

where Tq =8.5 kpc denotes the distance between the sun and the Galactic Center. The 
vertical distribution of SNR is taken to be box-shaped with a half thickness of 150 pc. 

Our model allows us to use true 3D source distributions. We know that the Galaxy 
has structure in the form of spiral arms, a bulge and so forth. We have thus folded a spiral 
arm model QGeorgelin and Georgelin 1976 , Vallee 1995| ) with the radial SNR distribution 



as described above to investigate the influence of spiral arm structure. In the Georgelin 
and Georgelin model the Galaxy has four symmetric arms with pitch angle of around 12°. 
We have rescaled their model to Tq =8.5 kpc. Each spiral arm is described by a Gaussian 
of dispersion 500 pc, i.e. a FWHM of 1177 pc. The spiral arm model is normalized in 
azimuth, so that the integral Jq^ r dcp of that model yields unity, and then folded with the 
radial SNR distribution. The normalization is required to preserve the radial distribution 
of SNR. A face-on view of the resultant source distribution is shown in Fig.|^. 



4. Results 

In this section we show results for two choices of the spatial distribution of sources, 
the pure SNR distribution and the SNR distribution folded with a model of the spiral arms 
in the Galaxy. We will also show results for two choices of injection index, at first a fixed 
index s=2.0 for all sources, and then a normal distribution of indices with mean 2.0 and 
dispersion 0.2. We will not vary any other parameter in the propagation model, and stick to 
the best fit values given by Webber, Lee and Gupta (1992). In a forthcoming paper we will 
extend our model to nucleons and determine the propagation parameters self-consistently 
which fit the 7-ray data and the local spectra of primary and secondary cosmic rays. In 
this paper our emphasis is to show that an injection index of around s=2.0 for cosmic ray 
electrons is sufficient to explain the observed spectrum of diffuse high-energy 7-rays , while 
leaving all other parameters unchanged. 
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4.1. The local electron spectra 

We have shown in Sec. 2 that if electrons are accelerated in SNR, their local spectra 
above ~ 30 GeV will strongly depend on time, so that the direct electron measurements 
above this energy may deviate from the average electron spectrum. As a result, we are not 
required to choose the electron injection index according to the directly observed electron 
spectrum above ~ 30 GeV. Below a few GeV, on the other hand, the locally observed 
spectra are strongly affected by solar modulation, for which we do not have reliable models, 
so that only over roughly one decade in energy does the local electron spectrum provide 
clear data. The data from direct electron measurements compared with the steady-state 
spectra of our model are shown in Fig.^ The model spectra fit the data reasonably well in 
the relevant energy range up to 30 GeV. 



4.2. The 7-ray emission 

As an example we show in Fig.^ the spectra of the bremsstrahlung and inverse Compton 
emission in the direction of the inner Galaxy for the case of sources with injection indices 
following a normal distribution of mean s=2.0 and dispersion fi=0.2. The figure includes 
the results of our model, the observed spectrum, and a template of the 7r°-decay spectrum 
( Permer 1986|) . At around 5 GeV the intensity due to leptonic processes is higher than that 



due to hadronic vr" decay. Our model assumes that the power-law behaviour of the electron 
injection spectra persists to ~20 TeV. The high-energy 7-ray spectrum will rather directly 
reflect structure in the electron source spectra. If for example the true source spectrum 
would deviate from a simple power-law above 1 TeV, the Inverse Compton spectrum would 
show corresponding features above 50 GeV. 

When we consider the latitude distribution of the 7-ray intensity above 1 GeV, which 
is done in Fig.^, we find a high level of agreement. The fraction of the total diffuse intensity, 
which is due to leptonic emission, is almost constant between b=0° and b=10°. It goes up 
from ~6% in the Hunter et al. model, ~10% in our model for an injection index of s=2.4, 
to 30% - 48% for an injection index of s=2.0, enough to explain all the 7-ray excess. 

We can compare the longitude distribution of the observed diffuse 7-ray emission above 
1 GeV to that of the leptonic contribution in our model. This is done in Fig.|^. It is obvious 
that the galactocentric gradient of the leptonic 7-ray emission in our model is stronger than 
in the data. This is the case for all cosmic ray propagation models which are based on the 
SNR distribution ([Webber, Lee and Gupta 1992| , pioemen et al. 1993| ). While towards the 



inner Galaxy the fraction of the total diffuse intensity, which is due to leptonic emission, is 
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around 35% - 52%, far enough to explain all the 7-ray excess, towards the outer Galaxy the 
leptonic contribution in our model accounts only for roughly two thirds of the excess. 

The SNR distribution of Case and Bhattacharya has a zero at galactocentric radius r=0 
which causes the double peak structure in the model intensity towards the inner Galaxy. 
Such double peak structure is not visible in the data, so that this effect may result more 
from the specific choice of mathematical function in the fit of the SNR distribution than 
astrophysical reality. We can generally say that a flatter SNR distribution would beget a 
better harmony of the longitude distribution of observed emission and model. Interestingly, 
one study of the SNR distribution, which was not based on a specific form of the radial 
profile, indicates very long radial scale lengths up to 9 kpc (|Li et al. 1991|). 




We have therefore investigated the impact of the fit uncertainties in the SNR 
distribution. Within 1 a in both parameters the SNR distribution may be 

-2.96 ^ I (21) 

r& J 

Here we show the results for the spiral arm model only. For similar local electron spectra the 
leptonic contribution to the diffuse Galactic 7-ray emission above 1 GeV can be 33% -47% 
while the center/anticenter contrast in the model agrees with the observed one to better 
than 25%. As can be seen in Fig.|, the latitude distribution of the leptonic contribution 
varies little. The longitude distribution of the leptonic contribution in Fig.^ shows that the 
overall gradient is reasonably well reproduced, but the double hump structure towards the 
inner Galaxy remains as well as a general overprediction around |Z| ~ 45°. This double 
hump structure however is a consequence of the mathematical function chosen to fit the 
SNR distribution and not a consequence of the SNR distribution itself. 



4.3. The synchrotron emission 

A further constraint on our model is the synchrotron flux towards the North Galactic 
Pole (NGP). Available data at 408 MHz ( [Haslam et al. 1982| ), at 820 MHz ( |Berkhuijsen 



1972| ), and at 1420 MHz ( [Reich and Reich 1986| ) can be corrected for zero level uncertainties. 



the contributions of the microwave background and unresolved extragalactic sources ( [Reich 



and Reich 1988[) , and contributions from the Coma cluster ( ^chlickeiser, Sievers and 



Thiemann 1986[) . Here we do not use data at frequencies below 100 MHz, which we expect 



to be affected by free-free absorption. In Fig.|T^ we compare the synchrotron intensity 
predicted by our model with the data in direction of the NGP. It can be seen that for a 
total magnetic field strength of lOyuG there is good agreement. 
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In our model the FWHM of the z-distribution of synchrotron emission at 1420 MHz is 
4.1 kpc, which is the value typically found for edge-on galaxies ( piummel 1991| ). 



5. Summary and discussion 

In this paper we have investigated whether cosmic ray electrons can be responsible for 
the recently observed high intensity of diffuse Galactic 7-ray emission above 1 GeV. Models 
based on the locally observed cosmic ray spectra underpredict the observed intensity by 
nearly 40% ( [Hunter et al. 19"^ ). One feature of these models is the relatively soft electron 
injection spectral index of s=2.4 ( |Skibo 1993| ), which is required to account for the local 
electron spectrum above 50 GeV. 

The recent detection of non-thermal X-ray synchrotron radiation from the four 
supernova remnants SN1006 ( [Koyama et al. 1995| ), RX J1713. 7-3946 ( p<.oyama et al. 1997| ), 



IC443 ( [Keohane et al. 1997| ), and Gas A ( [Allen et al. 1997| ) supports the hypothesis that 



Galactic cosmic ray electrons are accelerated predominantly in SNR. We have shown in this 
paper that, if this is indeed the case, the local electron spectra above 30 GeV are variable 
on time scales of a few hundred thousand years. This variability stems from the Poisson 
fluctuations in the number of SNR in the solar vicinity within a certain time period. While 
the electron spectra below 10 GeV are stable, the level of fluctuation increases with electron 
energy, and above 100 GeV the local electron flux is more or less unpredictable. 

With that time variability in mind we have seen that an electron injection index 
of s=2.0 is consistent with the data of the direct particle measurements if SNR are the 
dominant source of cosmic ray electrons. In fact, both the radio spectra of individual SNR 
( [Green 1995[) and the hard spectrum of the inverse Gompton emission at high latitudes 
( Phen, Dwyer and Kaaret 1996[) would better harmonize with an injection index of s=2.0 



instead of s=2.4. 

We have then presented a three-dimensional steady-state diffusion model for cosmic 
ray electrons, based on the propagation parameters which have been derived from similar 
models for cosmic ray nucleons. While being entirely consistent with the local electron 
electron flux up to ~ 30 GeV energy and with the radio synchrotron spectrum towards the 
North Galactic Pole, the leptonic contribution to the diffuse Galactic 7-ray emission above 
1 GeV in the Galactic plane increases from ~6% in the model of Hunter et al., ~, 10% in 
our model for an injection index s=2.4, to 30-48% for an injection index s=2.0 depending 
on the assumed spatial distribution of SNR and depending on whether some dispersion of 
injection spectral indices is allowed. An electron injection index s=2.0 can therefore explain 
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the bulk of the observed 7-ray excess over the predictions of the Hunter et al. modeL 

While the latitude distribution of the leptonic 7-ray emission is fully consistent with 
that of the observed emission, we find that the longitude distribution deviates from the 
observed one. In our model the contrast between the Galactic Center and the anticenter 
is stronger than in the data. A similar effect can be found in all cosmic ray propagation 
models which are based on the SNR distribution ( [Webber, Lee and Gupta 1992| , [Bloemeii 



et al. 19^ ). Note that in the Hunter et al. model only the emissivity spectrum is taken 



according to a propagation calculation ( [Skibo 1993|) , while the spatial distribution of the 



emissivity is scaled according to the distribution of thermal gas. 

We have seen that the fit uncertainties in the SNR distribution of Case and 
Batthacharya allow us to use a flatter profile, which leads to a better agreement between 
model and data in the longitude distribution (see Fig.^. Thus the gradient problem may 
be simply the result of an inappropriate choice of radial SNR distribution or lack of error 
propagation, respectively. This fiatter profile would also harmonize better with the results 
of Li et al. (1991), who found very large radial scale lengths of up to 9 kpc for the galactic 
distribution of SNR. Other possible sources of systematic errors are discussed below. 

This gradient problem is unlikely to be caused by additional thermal matter. Very 
cold (3°K) molecular gas has recently been discussed as candidate for baryonic dark matter 
( Pfenniger et al. 1994|) . If organised in small clumps (|Pfenniger and Combes 1994]) the 
probability of finding absorption features in the spectra of bright background objects 
would be small and the clumps would easily evade detection (for a review see [Combes and| 
[Pfenniger 1997] ). However, the thermal gas mainly affects the bremsstrahlung and only 



indirectly the inverse Compton emission, which dominates above 1 GeV 7-ray energy, and 
thus has little influence on the gradient in the total leptonic 7-ray emission above 1 GeV. 

If our model for the interstellar radiation fleld were wrong, then it would have only 
a limited influence on the gradient, since the Galaxy acts as a fractional calorimeter at 
high electron energies ( [Pohl 199^ , which channels the electron source power directly into 
synchrotron and inverse Compton emission with a ratio corresponding to that of the energy 
densities in the magnetic fleld and the photon fleld. Any radial variation of the magnetic 
fleld strength will be directly reflected in the Center/anticenter contrast of the synchrotron 
intensity, so that the radio surveys would constrain the parameter space here. 

It is our personal view that the uncertainties in the radial distribution of SNR are 
large enough so that the overly strong gradient in our model may simply be the result of an 
inappropriate choice of SNR distribution. The artiflcial double peak structure towards the 
inner Galaxy is an example of systematic effects which arise from possibly ill-deflned flts to 
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the SNR distribution. 

Our findings indicate a potential problem in the determination of the extragalactic 
7-ray background (preekumar et al. 1998|) . The standard method uses a linear regression 
analysis of observed intensities and model predictions (with the model of Hunter et al.) 
to extrapolate to zero Galactic intensity. If at higher 7-ray energies the inverse Compton 
emission is indeed much stronger than assumed by Hunter et al., then a large fraction of it 
will be attributed to the extragalactic 7-ray emission. Our steady-state model predicts an 
intensity of high-latitude inverse Compton emission above 1 GeV at a level of ~ 40% of the 
extragalactic background intensity, while Hunter et al. assume much smaller values. The 
intensity difference between the inner Galaxy and the outer Galaxy at medium latitudes 
(40° — 50°) would be around 10% of the extragalactic background intensity, depending on 
the choice of model for the electron source distribution (SNR distribution). On the other 
hand, if we are indeed living in a region of temporarily low flux of high-energy electrons, we 
would also expect the intensity of inverse Compton emission towards the Galactic Poles to 
be less that the steady-state value, so that the true level of Galactic 7-ray intensity at high 
latitudes is difficult to assess. 

The spectrum of inverse Compton emission, s~1.85 in our model, is somewhat harder 
than the s~2.1 of the extragalactic background. It has been noted before that the average 
7-ray spectrum of the identified AGN is softer than that of the background (Pohl et| 
al. 1997b| ), which indicates a problem with the idea that the background is mainly due 



to unresolved AGN. Now, since the presently determined background spectrum may be 
substantially contaminated with hard inverse Compton emission, the true background 
spectrum may be softer than s=2.1, which would probably relieve the spectral discrepancy 
with the average spectrum of resolved AGN (|Pohl et al. 1997b|). Anyway, the systematic 



uncertainty in the spectral index of the extragalactic 7-ray background emission is much 
higher than the statistical uncertainty. 

In a forthcoming paper we shall discuss the distribution of synchrotron emission in 
more detail. We shall also use a truly three-dimensional calculation of the interstellar 
photon fields based on COBE/FIRAS data. Finally we shall discuss cosmic ray nucleons in 
parallel with electrons to derive the propagation parameters self-consistently. 
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Fig. 1. — The locally observed electron spectra in the upper panel compared with the range 
of possible spectra in our model. The parameters of the model are given in the main text. 
For each experiment the 1 a uncertainty range is indicated by a grey-shaded band which 
connects the data points at the mean energies of the corresponding energy bins. The scatter 
between the results of different experiments indicates the level of systematic uncertainties. 
The range of possible spectra in our model is given by the grey-shaded bands in the lower 
panel. During 68% of the time the locally observed spectra will be in the dark grey shaded 
region, and during 95% of the time they will be within the light grey shaded region. The 
dashed black curve gives the 68% range for a weaker energy dependence of the diffusion 
coefficient (a=0.33 instead of a=0.6) to show the influence of this parameter. The white 
dash-dotted line shows one of the 400 random spectra as a particular example of what may 
be observed. The white dotted line indicates the time-averaged spectrum. The effect of 
solar modulation is taken into account for all model spectra using a force-field parameter 
$ = 400 MV ( pleeson and Axford 19681) . The data are not in conflict with the range of 
possible spectra in our model. 



-27- 




Fig. 2. — Face-on view of the surface density of gas in the Galaxy. The square indicates the 
position of the sun. The color scale is linear between surface mass densities of -2 MqPc~'^ 
and 30 MqPC~'^. The plot includes atomic, molecular, and ionized gas integrated from -500 
pc to 500 pc height above the plane. 



-28- 




Fig. 3. — Face-on view of the SNR distribution in the Galaxy according to the spiral arm 
model. The square indicates the position of the sun. The color scale is linear between zero 
and peak value. Note that the sun is located in an interarm region. 



- 29 - 




Fig. 4. — The data for the local flux of cosmic ray electrons compared with the steady 
state model spectrum. Solar modulation is taken into account using a force-field parameter 
$ = 400 MV (IGleeson and Axford 1968|) . For the case of the pure SNR distribution as 
electron source distribution, the solid line is for sources with fixed injection index s=2.0, while 
the dotted line is for source indices with a normal distribution of mean 2.0 and dispersion 
0.2. When the SNR distribution in spiral arms is taken as source distribution, we obtain the 
dashed line for sources with fixed injection index s=2.0 and the dash-dotted line for source 
indices with a normal distribution of mean 2.0 and dispersion 0.2. 
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Fig. 5. — The 7-ray intensity in the direction of the inner Galaxy. The data points are 
taken from Hunter et al. (1997). The error bars include an estimate for the systematic error 
of 8%, which accounts for the uncertainty in the energy- dependent correction of the spark 
chamber efficiency ([Esposito et al. 1998| ). The data are compared with bremsstrahlung 
(ebr) and Inverse Compton (ic) spectra from our model, here on the basis of sources with 
injection indices following a normal distribution of mean 2.0 and dispersion 0.2 and the 
spatial distribution of SNR in spiral arms. The vr'^-component is a template and not a 
model. 
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Fig. 6. — The latitude distribution of the 7-ray emission above 1 GeV compared with the 
model prediction for the leptonic contribution. Again the solid line is for fixed injection 
index, and the dotted hne is for a distribution of indices with dispersion 0.2. The numbers 
give the percentage of the observed emission in certain directions which is due to leptonic 
processes. The top panel displays results for the pure SNR distribution as electron source 
distribution, while the bottom panel is for a source distribution of SNR in spiral arms. 
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Fig. 7. — The longitude distribution of the 7-ray emission above 1 GeV. The upper soUd 
Hne shows the observed distribution according to Hunter et ah (1997). The distribution of 
the leptonic emission in our model is shown for comparison, here for the SNR distribution 
as source distribution. Again the solid line is for fixed injection index, and the dotted line 
is for a distribution of indices with dispersion 0.2. The numbers give the percentage of the 
observed emission in certain directions which is due to leptonic processes. The top panel is 
for the pure SNR distribution and the bottom panel is for a source distribution of SNR in 
spiral arms. 
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Fig. 8. — The latitude distribution of the 7-ray emission above 1 GeV as in Fig.^ except that 
here the flatter SNR distribution in spiral arms has been used. The indicated percentage of 
leptonic contribution of 33% - 47% is sufficient to account for all the excess. 
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Fig. 9. — The longitude distribution of the 7-ray emission above 1 GeV as in Fig.|^ except that 
here the flatter SNR distribution in spiral arms has been used. The overall gradient is well 
reproduced. The double hump structure towards the inner Galaxy and the overprediction 
around |Z| ~ 45° indicate that the mathematical profile in the fit to the SNR distribution 
may be ill- defined. 
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Fig. 10. — The synchrotron intensity in the direction of the North Galactic Pole. The error 
bars indicate data from surveys with reasonable zero level cahbration. The solid line is the 
model spectrum based on the pure SNR source distribution while the dotted hne is for the 
SNR distributed in spiral arms. 



